Influence of dissipation on the extraction of quantum states via repeated 

measurements 
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A quantum system put in interaction with another one that is repeatedly measured is subject to 
a non-unitary dynamics, through which it is possible to extract subspaces. This key idea has been 
exploited to propose schemes aimed at the generation of pure quantum states (purification). All 
such schemes have so far been considered in the ideal situations of isolated systems. In this paper, we 
analyze the infiuence of non-negligible interactions with environment during the extraction process, 
with the scope of investigating the possibility of purifying the state of a system in spite of the 
sources of dissipation. A general framework is presented and a paradigmatic example consisting of 
two interacting spins immersed in a bosonic bath is studied. The effectiveness of the purification 
scheme is discussed in terms of purity for different values of the relevant parameters and in connection 
with the bath temperature. 

PACS numbers: 03.65.Ta, 03.65.Yz, 32.80.Pj, 32.80.Qk 



I. INTRODUCTION 

During the last decades, enormous progresses have 
been made in such physical contexts as cavity quantum 
electrodynamics (CQED) m, superconductor-based cir- 
cuits [1], and trapped ions Q. In connection with the ap- 
plications in the field of quantum information Q , seminal 
experimental results have been reached, such as the im- 
plementations of quantum logical gates 0, 0] and the 
realizations of quantum teleportation Generation of 
quantum states is a crucial issue in nano-technologies, be- 
cause it is the basis of initialization processes with which 
many experimental protocols begin. Incidentally, it also 
provides the possibility of observing physical systems be- 
having according to the predictions of quantum mechan- 
ics once a nonclassical state has been created. 

Recently, a strategy for the generation of pure quan- 
tum states through extraction from arbitrary initial 
states has been proposed d, [l^l ■ This procedure is based 
on the idea of putting a quantum system in interaction 
with another one that is repeatedly measured in order to 
induce a non-unitary evolution which forces the former 
system onto a Hilbert subspace. If such a subspace is 
one-dimcnsional, the process reduces to the extraction of 
a pure quantum state. For this reason, this procedure has 
been addressed as a "purification" [j, [ll| . On the basis 
of this general scheme, many applications have been pro- 
posed: it is possible to extract entanglement ^^'^ 
the initialization of multiple qubits would be useful for 
quantum computation |lOl . [l3| ; extensions of the scheme 
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enable us to establish entanglement between two spa- 
tially separated systems via repeated measurements on 
an entanglement mediator [l^; in single trapped ions, the 
extraction of angular-momentum Schrodinger-cat states 
has been proposed and the possibility of steering the 
extraction of pure states through quantum Zeno effect 
has been predicted (l6j . In passing, we mention that 
the approach here recalled is related to quantum non- 
demolition measurements [l7j . introduced for gravita- 
tional wave detection and exploited in different physi- 
cal systems for applications in quantum computation and 
information [l^ and for quantum state generation in gen- 
eral. For instance, in the context of trapped ions, putting 
the vibrational degrees of freedom in interaction with the 
electronic degrees of freedom, and repeatedly measuring 
the atomic state, it is possible to generate Fock states 
both in one-dimensional [20| and two-dimensional con- 
texts [m. 

Until now, the extraction of pure states through re- 
peated measurements has been considered in ideal situa- 
tions, that is, the evolution is assumed to be unitary ex- 
cept for the change of state by the measurements. In this 
paper, we discuss how the predictions change when the 
system is put in interaction with an environment and, as 
a consequence, is subject to a non-unitary evolution (even 
between two successive measurements) , which is assumed 
to be described by a Lindblad-type equation [23|. First 
of all in this paper, we discuss the behavior of a quan- 
tum system, whose dynamics is governed by the repeated 
measurements represented by projection operators on the 
other system in interaction with the former and by the 
dissipative environmental interaction with the two quan- 
tum systems. A criterion for the extraction of pure states 
in the presence of dissipation is derived and a quantum 
system composed of two two-level systems (qubits) im- 
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mersed in a common bosonic bath is analyzed as a very 
simple example. The ability or disability of extracting 
pure states for a one qubit system by this scheme is esti- 
mated numerically in terms of the purity of the density 
operator. 

The present paper is organized as follows. In Sec. [Tll 
the general aspects of the purification protocol in the 
presence of interaction with an environment arc dis- 
cussed. In Sec. mil a simple model of two mutually in- 
teracting two-level systems immersed in a bosonic bath 
is analyzed, and the main differences between the ideal 
situation (in the absence of dissipation) and the more 
realistic situation here studied (in the presence of dissi- 
pation) are singled out. Finally, some conclusive remarks 
are given in Sec. IIVI 



II. GENERAL FRAMEWORK 



Framework 



Assuming that the time elapsed between two successive 
measurements is always r and that every of the A'^ -I- 1 
measurements has confirmed X to be in the state |$)x, 
the initial state /Ctot(O) is mapped into the final one. 



ptotir, N)^VM{t) 
where 



(3) 



Vtotir) =VM{t)V, 



(4) 



which, although rigorously speaking acts on Liouville 
space of the whole system S-l-X, is substantially a 
superoperator acting on the Liouville space of subsys- 
tem S, except for the trivial action on X (i.e. projection 
on the state |<I')x)- Therefore, it can be rewritten as 
Vtot(''') = V{t)(E)'P, where V(t) is a superoperator acting 
only on the Liouville space of S. 

If the master equation is given in the Lindblad 



form [23, the evolution of S-l-X in interaction with the 
bath is written down in the Kraus representation [23] as 



The scheme of extraction wc study in this paper is 
based on the idea presented in Rcf. Q . Assume that we 
are interested in preparing a target system, S, into a pure 
state. To this end, put it in interaction with an ancilla 
system, X, which is repeatedly measured and found in 
the same state, say |$)x- As a consequence of both the 
interaction with and the measurements on X, system S 
is subject to a non-unitary dynamics which forces it in a 
subspace. If the subspace is one-dimensional, the relevant 
state is extracted. We underline that if a negative result 
is obtained measuring system X (i.e., if the ancilla system 
is found in a state different from the expected one) , then 
the trial can be considered as "failed" and the system 
has to be reset in order to restart the experiment, or 
we can just continue the process as a new trial with the 
initial condition given by the result of the unsuccessful 
projection. 

In the following, we consider this scheme but in the 
presence of dissipation, in order to examine whether a 
pure state can still be extracted. To this end, assume now 
that both S and X are immersed in a bath, so that the 
dynamics of the whole system S-t-X between the repeated 
measurements is governed by the master equation 



dpt. 



dt 



Cpt. 



C = H + V, 



(1) 



where ptot is the density operator of S+X (the trace over 
the bath degrees of freedom has already been performed) . 
Ti. is the Hamiltonian superoperator defined by Tiptnt = 
— Ptot], while V is the Lindblad-type dissipator [23 
which takes into account the interaction with the environ- 
ment, in the Markovian limit. Let M.{t) = e'~'^ denote 
the superoperator describing the dissipative dynamics of 
S-l-X and V the projection superoperator associated with 
the measurement of the state |$)x, i-C 



Vptot = |$)x($|ptot|$)x($|. 



M{T)ptot =Y.n{T)ptotfl{T), 



(5) 



where G -I- 1 is the number of Kraus operators Tfc (t) in- 
volved in the decay process. Therefore, the superoperator 
V(t) that maps the state of S just after a measurement 
to another after the next measurement reads 



V(r)p- Vl4(r)pl4t(T), 



fc=0 



where p is the state of S and 

Vfe(T) = x($|rfc(T)|$)x 



(6) 



(7) 



is an operator acting on the Hilbert space of S. Note that 
the state of the whole system S-l-X after a measurement 
is factorized like ptot = p® |$)x('I*|- 

The properties of V(r) determine the fate of the state 
of S, after many repetitions of measurements on X. In 
fact, assume that V(t) is diagonalizablc and consider 
its spectral decomposition in terms of its eigenprojec- 
tions IH: 



V(r)-^A, 



(8) 



where n„ is the eigenprojection belonging to the eigen- 
value A„ and satisfies the ortho-normality and complete- 
ness conditions 



Ilmlln — ^7nn^^m 



1. 



(9) 



The eigenvalues A„ are complex-valued in general and 



ordered in such a way that |Ao 



|Ai| 



I As 



> 



> lA 



so -1-2 1 



> 



If V(t) is not diagonalizable. 



(2) the Jordan decomposition applies instead of 
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Generalization of the following argument to such cases 
is st raig htforward. See, for instance, an Appendix of 
Ref. [ifll. 

The evaluation of the A^th power of V(r) shows that, 
the larger the number N of the repeated measurements is, 
the more dominant the blocks belonging to the maximum 
(in modulus) eigenvalues are over the other blocks. That 
is, for a large enough N, the action of V^(t) diminishes 
the components of the system density operator which do 
not belong to the generalized eigenspaces corresponding 
to the maximum (in modulus) eigenvalues. Indeed, 



V^(r)=5:A,^n„ 



90 

E 

T1=0 



(10) 



which, in the case wherein only one eigenvalue exists 
whose modulus is maximum (i.e. when go ^ 0), reduces 
to 



N- 



(11) 



Therefore, in such a case, for a large enough (depending 
on the structure of the spectrum) N, the action of V^(r) 
essentially reduces to the action of Hq, hence projecting 
the system in the relevant subspace. 

As for the outcome state, it is worth mentioning that, if 
Ao is degenerated, the state of S extracted by XIq depends 
on the initial state of S, p(0), since such a final state is 
substantially proportional to nop(O), which depends on 
p(0) if XIq refers to a multi-dimensional space. On the 
contrary, if Aq is not degenerated and Hq refers to a one- 
dimensional subspace, i.e. 



nop = fo{p)PQ 



(12) 



with pq the relevant eigenstate (which is an element of 
the Liouville space of S) and /o a suitable form on the 
Liouville space of S, the eigenstate po is extracted irre- 
spectively of the initial condition of the system, provided 
nop(O) 7^ 0. As a first step, in this paper we shall con- 
centrate on such situations wherein there is a single ex- 
tracted eigenspacc which in addition is one-dimensional. 

In an ideal case wherein any sources of dissipation are 
absent, Eq. ^ reduces to V{t)p = Vo{t)pVq{t) with 

Vq{t) = x('I'|e"*^*|$)x 3, so that, if Vo{t) possesses a 
nondegenerate and unique maximum (in modulus) eigen- 
value, then the final state to be extracted is the pure 
state |wo)s(ito|, since in this case Eq. ([T2|) is supplemented 
by Po = \uo)s{ua\ and fo{p) = s{vo\p\vo)s where |uo) 
and {vo\ are the right- and left-eigenvectors of Vo{t) be- 
longing to its largest (in modulus) eigenvalue, respec- 
tively. This is the basic idea of the purification scheme 
based on the repeated measurements, which was first pro- 
posed in Ref. [9| and has been analyzed and developed in 
Refs. [l3,[il,[ll[i3,[ll,[ii. On the other hand, it is im- 
portant to stress that in the non-ideal case, even if there 
is a single and nondegenerate eigenvalue that is maxi- 
mum in modulus, this does not guarantee that the state 
Po is a pure state and hence there is no warranty that 



the final extracted state is pure, either. In this sense, our 
state-extraction scheme may not necessarily be an effec- 
tive purification scheme. However, we still try to seek 
a possibility of extracting a pure state even in the pres- 
ence of dissipation. The examination of such situations 
wherein we can extract a pure state is the main topic of 
this paper. 

Efficiency — This scheme for extracting quantum 
states is a conditional one, in the sense that each time 
system X is measured it has to be found in the same 
state, denoted by |$)x. In Refs. 0, [l^, [H, d, Q [H, [3 
it is proved that the probability of success of the ex- 
traction, that is the probability of finding system X 
in the state |$)x successively A^ times, is given by 
the normalization factor of the state extracted by ([3]), 
trs{V^(r)p(0)} = E„A^trg{n„p(0)}, which behaves 
asymptotically as — > A^ trs{nop(0)} as A^ cx) [or 
— > J2n<go trg{n„p(0)} in the more general situation, 
which from now on we shall not mention anymore for the 
sake of simplicity] . These expressions for the probability 
of success (still valid in the non-ideal case, provided the 
projectors n„ are the appropriate ones) show that the 
structure of the spectrum of V(t) plays a crucial role for 
the efficiency and fastness of the extraction. In partic- 
ular, on the one hand, the fact that Aq has a modulus 
quite larger than those of the other eigenvalues makes 
V^(t) quickly approach A^IIq, so that a smaller number 
of measurements is required to well approximate the fi- 
nal result llop{0). On the other hand, the closer to unity 
the modulus of Aq is, the greater the probability of suc- 
cess is, approaching just trs{no/9(0)} (without decaying 
out completely) for Aq — 1. On the contrary, for small 
values of Aq , the probability quickly approaches zero like 
A^ , which means that if a large number of measurements 
is required to approach the state (trs{no/c(0)})~^no/o(0) 
the scheme becomes very inefficient. Therefore, the num- 
ber of measurements necessary to extract the target state 
would be an important measure of efficiency. It can be 
roughly estimated by the following argument. The idea is 
to sec how much the relevant part approaching the target 
state, A^no/o(0), dominates over the rest: 



EJA„nin„p(o)|| 



(13) 



where 



is a certain norm, for instance \\A\\ 



■\/trg{At^}. This measure approaches unity in the limit 
of an infinite number of measurements, p 1 as N —^ oo. 
Then, we ask how many measurements are necessary for 
this quantity to exceed a desired value < po < 1- After 
rewriting m as A^||nop(0)||/(E„^o |A"ri|n"P(0)!l) = 
p/{l — p), a sufficient condition for p > Pq is given by 



Ao^||nop(0)|| 



> 



Po 



(M-l)|Airi?(p(0)) - 1-po' 



(14) 



where M is the number of eigenvalues or equivalently 
the dimension of the Liouville space and R{p{0)) = 



4 



maxn^o l|n„p(0)||. The number of measurements neces- 
sary to get a better quality than pq is therefore estimated 
by 



hibo/(l-Po)]+ln(A/ 

N > 

ln|Ao/Ai| 



1) 



, J?.(p(0)) 

'^^^ l|nop(o)|| 



(15) 



It is important to note that this threshold depends on 
p{Q), according to the expectation that the larger is the 
norm of the relevant part in the initial state, |jno/o(0)||, 
the smaller is the number of necessary measurements. 



B. Searching for Pure Eigenvectors: the Criterion 

It is easy to show that a necessary and sufficient con- 
dition for a pure state being an eigenvector of the map 
V(r) is that it is a simultaneous eigenstate of all the op- 
erators Vkir) involved in the relevant map (see Eq. ©). 
The proof of this statement proceeds as follows. 

=J> Obviously, if the state \(j))s is a common eigen- 
state of all Vfc's, i.e., Vl-|(/))s = o:k\4>)s and consequently 



then one has 



G 



Vkirmsm^ir) = XMsi^l (16) 



fc=0 



where 



= ^ a^ak > 0. 



(17) 



<j= Let the pure state \cj))s{4'\ be an eigenvector of 
V(t), then EtoMrMsmiir) = A^|</.)s(0|. Con- 
sider now the overlap with a quantum state \(j>±)s{4>±\ 
orthogonal to \<j>)s{4>\- 



o = A^s(0±i</': 

G 

k=0 
G 

k=0 



(18) 



from which it follows that s{4'±\^k{T)\(f>)s = for all k 
and whatever the state |c/>j_)s is, provided it is orthogonal 
to \4>)s- In other words, it means that 



H(r)|0)s =«fe|0)s, Vfc, 



(19) 



where ak is a suitable complex number. This completes 
the proof. 



C. Searching for Pure Eigenvectors: Purity 

In those cases in which we are not able to extract an 
exactly pure state, there is a possibility of extracting "al- 
most pure" states, that is, mixed states very close (in the 



sense of purity) to the pure states. To look for almost 
pure states which can be extracted, let us recall a measure 
of purity of a given state. We show later how the purity 
of the eigenstate of the linear map V(r) corresponding 
to the maximum eigenvalue behaves as a function of the 
parameters of the scheme, i.e. the interval of time r and 
the repeatedly measured state of X, |<I')x- 

The purity of a state is defined as the trace of the 
square of the relevant (normalized) density operator [26| : 



P{p) = trs 



(20) 



This quantity is upper and lower bounded in accordance 
with 1/L < P{p) < 1 with L being the number of levels of 
the system under scrutiny. Observe that the maximum 
value [P{p) = 1] corresponds to pure states, while the 
minimum value [P{p) = ^/L] corresponds to maximally 
mixed states with maximal von Neumann's entropy. 



D. Weak-Damping Case 

It is possible to derive a formula for the purity of the 
extracted state for general systems in the weak-damping 
regime. Such a formula would be useful for understand- 
ing which parameters spoil the purity and convenient for 
an optimization of the purification. 

Let us decompose the relevant map V(r) into two 
parts. 



(21) 



where V'''-'(r) is the map in the absence of the environ- 
mental perturbation and the rest is treated as a pertur- 
bation to it, which is given in the weak-damping regime 

by 

5V{t)®V~[ dtVe^^^-*^Ve^*V. (22) 
Jo 

Assuming that Vo(t) is diagonalizablc, let |u„)s and 
s{vn\ denote its right- and left-eigenvectors, respectively, 
which form a complete ortho-normal set, ^„ \un)s{vn\ = 
Is [131 • (We also normalize the right-eigenvectors as 
s{un\un)s = !•) Then, the right-eigenvectors of the ideal 
map read 

VW(r)a(°) =AWa(°), =|«,„)s(«„|. (23) 

These are orthogonal to the left-eigenvectors 

^^^l = \v,n)s{Vn\ (24) 

in the sense 

^Ll') = S„,m'5nn', {A, B) = tvs{A^B}. (25) 

We are interested in a situation where we can purify S in 
the absence of the environmental perturbation. That is, 
Aqq^ is not degenerated and is the only eigenvalue that is 
the largest in modulus. 



5 



Now, the standard perturbative treatment yields the 
first-order correction to the right-eigenvector, 



JO) {a^l„SV{r)a^nl) 



7n' n' ^7nn 



.(0) 



(0) 



(26) 



where the constant c is set equal to zero by the normal- 



ization condition (tTmn + Samn, (Tmn + Samli) = 1. 



' mn 
i(0) 



This formula is valid when Amn is not degenerated. 
We are interested in the state to be extracted, i.e. po — 



00 

)/(l + trs Saoo 



(1) 



Since Aqq^ has been assumed 
(1) 



to be nondegenerated, the formula ([26|) is valid for JcTpo 
The purity of po up to this order is therefore given by 



PiPo) 



1 



2 (trs Sa^Qo' - s(ito|<5(T^oVo>s) 



^m/ S 

mn/00 



^00 



i(0) 



(27) 



where Qo = Is — |'"o)s(iio| is a projection operator. 

This is the formula for the purity of the extracted state 
up to the first order in the decay constants in the weak- 
damping regime. This shows that, if the state CTqq^ to be 
extracted in the ideal case is an eigenstate of the pertur- 
bation, i.e. 5V(T)crQQ^ cx o-QO^ ll''^ first-order correction to 
the purity vanishes and the purification is robust against 
the environmental perturbation, at least up to this or- 
der. This is a weaker version of the criterion discussed in 
Sec. Ill Bl and is convenient since the dissipator of a master 
equation, 2?, suffices to this criterion without knowing the 
Kraus operators Tfc(r) of the decay process, which may 
require solving the master equation. Furthermore, this 
formula would be useful for finding a parameter set that 
optimizes the purity (minimizes the first-order correction 
to the purity). 

When S is a two-level system, the formula (P7)) is re- 
duced to 



P(po) ^ 1 - 



x(0) 



(28) 



III. A SIMPLE MODEL 

In this section, we apply the ideas presented above to 
the case of a simple model. Such a system consists of two 
mutually-interacting spins immersed in a bosonic bath, 
one of which is repeatedly measured to purify the other. 



A. Model 

Two-spin system — Consider a system of two in- 
teracting spins or pseudo-spins, for instance a couple of 



identical two-level atoms subjected to a dipolar coupling. 
Assuming that the matrix elements of the dipole opera- 
tors are real, and neglecting the counter-rotating terms, 
one reaches the following Hamiltonian (for details, see 
Refs. EH!): 



^ 2" 

i=S,X 



(1 



(S)-(X) 



JS)-(X) 



) (29) 



where a« = |T).(T| - U).(i|, ^V'' = |T).(il = i^'-'V, ^ 
is the Bohr frequency of the two-level system and e the 
coupling constant. We have set ?i = 1. 

The eigenstates of the Hamiltonian are the triplet and 
singlet two-spin states: 



|2)tot = |T)slT>x, 

|l)tot = ^[|T)sU)x + U)s|T)x] 

|0)tot = U)sU)x, 

1 



^LIT)sU)x 



U)slTb 



(30a) 
(30b) 
(30c) 
(30d) 



which are common eigenstates of and with S = 
^(o-'-^' 4- (T*^^^), whose eigenvalues are given by E(S -I- 1) 
and TOs, respectively. The corresponding eigenenergies 
are 2il, ft + e, 0, and fJ — e, respectively. If we consider 
the case > e, then |0) is the ground state. 

Interaction with a bosonic bath — The interaction 
with a bosonic bath, whose free Hamiltonian is given by 
Hb = ^ dkujka\ak, is modelled through the system-bath 
interaction Hamiltonian 



E 

i=S,X 



(i?,+i?|)(a«+a«). 



(31) 



where Bi = J dk gk{ri)ak, with the position of spin 
i and gki^i) the coupling constant between the atom at 
position ri and bath mode k. Following the standard 
derivation [2^ and assuming the spins very close each 
other in order to have 5fc(ri) ~ gk{f2): we reach the 
following master equation in the Schrodinger picture for 
the density operator ptot of S-l-X: 



dpto^ 
dt 



-i[iftot,Ptot] 

+ 72(1 + 71_)I?21/0tot -I- 7l(l + n+)l?ioptot 

+ 72"-2?i2Ptot + 7in+X'oiPtot (32) 



with Vijptot 



|j)tot(i|ptotK)tot(j| - \{\i)tot{i\:Ptot}, n± 
the mean numbers of bosons in the bath modes of fre- 
quencies ri ± e, which are the Bohr frequencies between 
the states involved in the transitions |2)tot |l)tot 
[Vt - e) and |l)tot ^ |0)tot (P- + e)- 7i and 72 are 
the decay rates related to such modes evaluated as the 
spectral correlation functions of Bi + Bj, and are re- 



2n J dk |gfep(5(cjfe — + e) and 
Finally, Htot is the 
Lamb-shifted Hamiltonian of S-l-X. 



lated to gfc's by 71 
72 = 27r / dk \gk\'^S{uJk - fl - e 
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B. Extraction of Pure States under the Influence 
of a Zero- Temperature Bosonic Bath 

Consider now the special case wherein the bath is at 
zero temperature. The spin labehed with X is repeat- 
edly measured and found in the state |<&)x = cos 5|t)x + 
e''<^ sin |||)x, while the other spin, labelled with S, is 
driven toward a quantum state through its interaction 
with X. The same situation is discussed in Ref. in 
the absence of the environmental coupling, where it is 
found that the extracted state can be made pure very 
efficiently, in particular measuring the states ||)x and 
||)x- The evolution of the damped system between two 
successive measurements is easily evaluated, for instance, 
following the approach developed in Ref. [30[ , and is given 
by (see Eq. [5]) 



ptot(t) = Y,Tk{t)ptomTi{t) 



(33a) 



fc=0 



with 4 Kraus operators. 



%{t) = |0)tot(0| +e-**e-'("+^)*|l)tot(l| 

(33b) 

fi(i) - Vl-e-Ti*|0)tot(l|, (33c) 



T2{t) = 



72 



71 - 72 



.(e-72t_e-7it)|l),„t(2|, (33d) 



m) = Jl + |0)tot(2|. (33e) 

V 7i ~ 72 



To(t) reduces to the unitary evolution operator in the 
case 72 = 7i = 0, whereas the others, i.e. Tk{t) for fc > 1, 
vanish. 

Now measure system X repeatedly every after t dur- 
ing the dissipative dynamics ([55)1 . According to (fTO|) . in 
order for the S state \(j))s{4>\, with \(j))s = cos^|t)s + 
e'^ sin §|i,)s, be a pure cigcnstatc of the contracted map 
V(t), it should satisfy 

s,(r) = s(^±|x($|rfe(t)|<i>)x|0)s =0, A; = 0,1, 2, 3 

(34) 

with \(j)±)s = sin ^|t)s — e'^ cos ||J,)s. Indeed, it is equiv- 
alent to look for the eigenstates of the contracted oper- 
ators Vkir) defined in ([7|). It is straightforward to find 
that 



Si (t) cx sin I cos ^ (e*^ cos | sin | + e*-^ sin | cos ^) , 
S2{t) cx cos I cos I (^e~^^ sin | sin | — e^'^ cos | cos |) 
S3 (r) (X sin | cos | sin | cos | , 



where the proportionality factors are the nonvanishing 
coefficients in p3cp - p3ep . From these expressions, it fol- 
lows that si(t) = S2(r) = S3(r) = is accomplished only 
for cos § = 0. This condition is necessary and sufficient 
to make Sfc's vanish for fc = 1, 2, 3. In order to make so(t) 
vanish too, it is necessary to have sin | = or cos § = 0. 
In fact, the condition cos | = means |0)s = |i)s, and 
evaluating so{t) in such a special situation provides 



so(T)=s(T|x($|To(T)|<i>)xU)s 



g~»Xe-'(n-e)r Q Q 

■sm-cos-(l-e e ), 
2 2^ ' 



V2 



(36) 



or cos 2=0- 



which; for t 7^ 0, vanishes only if sin 

This analysis shows that in some special cases, that is, 
when the state of X is repeatedly measured and found in 
|*)x = |T)x (sinf = 0) or |$)x = U)x (cos f = 0), the 
contracted linear map V(t) has the S state |i)s(i| as a 
pure eigenstate. To reach the final conclusion about the 
possibility of extracting such a pure state, one needs to 
know whether the corresponding eigenvalue is the max- 
imum (in modulus) in the spectrum of the map. We 
shall therefore diagonalize the contracted map in the two 



cases, sm ; 



and cos f = 0. 



Representing the density operator of S as a 
four-dimensional vector, p = {p^^ 1 Pili P^li Pi]) = 
(s(T|p|T)s,s(i|pU>s,s(T|pU)s,s(i|p|T)s), the contracted 
linear map V(r) is substantially represented by a 4 x 4 
matrix. 

Repeatedly measuring \[) x ((^ = t^) — In the case 
where the X state | j)x is repeatedly measured, the corre- 
sponding linear map V(r) = V|(t) is represented by the 
following matrix: 



Vi(r) 



/ l/i(r)P \ 

i(l-e-''i^) 1 

/x(r) 

V fliT)J 



(37) 



with/x(r) = ie-*(^^-^)-(H- 
ues of this matrix are 



) . The eigenval- 



Ao = l, Ai=A;=/i(r), A3 = |/i(7 



(38) 



(35) 



The right-eigenvector corresponding to the maximum 
eigenvalue Aq = 1 is the pure state po = |i)s(i|- The 
larger the time r is, the smaller the other three eigen- 
values of the map are and the faster the extraction of 
Pq = U)s(i| is, in the sense that it requires a smaller 
number of steps. 

Repeatedly measuring |t)jc (9 = 0) — In the case 
wherein the X state |t)x is repeatedly measured, the map 
V(t) reduces to 
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-72T 



72 (e 



VT(r) = 



2(71-72) 















'VT(r) 




(39) 



with /|(r) = ie-*("+'^)^(l + e-^^^e*2er)^ rpj^jg matrix 
is easily and exactly diagonalized as long as e~'^'^'^ ^ 
|/|(t)P. There are two cases in the ordering of its eigen- 
values. Case I: if e^^'^ < 



Case II: if e"^^ > |/|(t)|, 



(40) 



Ao = e-^^^ Ai =A* =e^^VTM, A3 = I/tM^. 

(41) 

In case I (which surely occurs in the strong-damping 
regime 72T 00), a pure state p*^^' = |i)s(i| is extracted, 
while in case II, a mixed state is extracted. 



with 



1 



PTTlT)s(TI+mU)sai 



Pli = 1 - PTT 



^ 72(e-^^"-e-^i") 

2(7l-72)(e-^^--|/T(r)P)■ 



(42a) 

(42b) 
(42c) 



The latter is not in contradiction with the previous state- 
ment that one has a pure eigenstate for sin | = 0. In- 
deed, the state |i)s(i| is still an eigenstate of the map, 
but it does not correspond to the maximum eigenvalue 
anymore, and then it is not the state to be extracted. 
The purity of the state p^^^' in ([^^ is given by 



1 



(43) 



In the weak-damping case 71T, 72T <C 1 (and assuming 
siner ^ for simplicity), one has a ~ ^ singer 
first order in 71 t and 72 t, and hence 



P(p(")) ~ 1 - 2a ~ 1 



I2T 
siv? er 



(44) 



This formula, that alternatively can be directly derived 
using (|28|) . shows that in the weak-damping regime (i) 
the purity is linearly affected by 72, while (ii) it is not in- 
fiuenced by 71. Furthermore, (iii) the purity is optimized 
by taking a nontrivial time interval r ~ 0.377r/e. 

It is worth noting that in the weak damping limit we 
cannot extract a pure state, while in the strong damp- 
ing limit a pure state can be obtained, which is the op- 
posite one would expect. To understand this fact, con- 
sider first of all that S-|~X has two stable states, |0)tot(0| 



and |s)tot(s| according to and second that in the 

strongly dissipativc case the system has time to relax 
onto the equilibrium state which is a mixture of the two 
stable states, whose statistical weights are determined 
by the initial condition. Then, repeatedly measuring the 
state |t)x cuts the population of |0)tot in the mixture and 
leaves only x(T|s)tot(s|T)x oc |i)s(i|- 

The case e~~'^ = \ f-\[T)\ corresponds to a degenerate 
case and hence, as clarified in Sec. |lTl is not in the scope 
of this paper since it does not permit the extraction of 
a precise state irrespectively of the initial state of the 
system. 

Notice that the general case corresponding to measur- 
ing a generic state |^>)x can be discussed in the weak 
damping limit through the perturbation analysis. 



C. At Finite Temperature 

The analysis on the model has so far been focused on 
the zero-temperature case and showed that the only pure 
state that can be extracted at the zero temperature is |J,)s 
when the state of X is repeatedly measured and found in 
|T)x or |J,)x. The question of what happens in the case 
of non-zero temperature naturally arises. To answer this 
question, we resort to numerical calculations. 

The linear map V(t) depends on r, the measured state 
|$)x (which is individualized in the Bloch-sphere by the 
polar and azimuthal angles, and x, respectively), and 
in general the temperature of the environment, T. Given 
the map, the eigenvector associated with the maximum 
eigenvalue, and its purity are functions of all such quan- 
tities (r, 9, X and T). 

In Fig. [U the purity of the state to be extracted is 
shown as a function of the parameters r and with x = 
being fixed, in four different situations concerning the 
bath. For the numerical calculations, we have set f2/e = 
10, 72/e = 0.1, and 71/72 = 0.95. 

In Fig. [IJa), the purity of the state to be extracted 
in the ideal situation is plotted, that is, in the absence 
of interaction with the bath. According to the discus- 
sion in Sec. [Tll the purity in such a case is expected to 
be equal to 1 whatever the parameters r, 6 and x are. 
In the other three figures, the behavior of purity in the 
presence of interaction with the bath at different tem- 
peratures is shown. Figure [T](b) refers to the case of very 
low temperature, effectively zero, and shows that in some 
regions of the parameter space, there is a possibility of 
extracting pure or almost pure states. This is not in con- 
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FIG. 1: The purity of the extracted state vs the parameters er and Qj-K (x = for all cases), in different situations: (a) 
in the ideal case, i.e. in the absence of interaction with the bath, and in the presence of interaction with the bath with (b) 
kBT/hO. — 0.01, (c) k^T/hO, = 1, and (d) ksT /liQ. = 10. In all cases, ratios between salient physical quantities are fixed as 
n/e = 10, 72/e = 0.1, 71/72 = 0.95. 



tradiction with the analysis in Sec. IIII Bl where it has 
been found that a pure state can be extracted only for 
= 0, TT. This result refers to an exactly pure eigenstate, 
while the numerical calculations here reported show the 
value of purity, which can be very close to unity although 
not exactly 1. 

In Fig.[T]Jc), one can see that in practice there is no re- 
gion in the parameter space corresponding to pure states: 
the purity is visibly smaller than unity everywhere. Fi- 
nally, in Fig. md), we see that at a higher temperature 
{k^T = IQMl, with fee the Boltzmann constant), the pu- 
rity of the state to be extracted is equal to the minimum 
value for the two-level system, i, almost everywhere, that 
is, irrespectively of the values of parameters t and 9. 

In Fig. [21 the purity is plotted as a function of the tem- 
perature T and of the time interval r between successive 
measurements, when a fixed state of system X character- 
ized hy 9 = 37r/4 and x = is repeatedly measured. It 
is well visible that the more the temperature increases, 
the more the purity of the extracted state approaches the 
minimum value, that is, \. 

All these results express in a clear way that the inter- 



action with an environment deteriorates the reliability of 
the purification scheme based on repeated measurements, 
although at the zero temperature pure states can still be 
extracted. 



D. Efficiency 

In a realistic situation, the probability of extracting the 
target state as well as the number of measurements one 
has to perform are important factors to consider. Accord- 
ing to the discussion at the end of Sec. Ill A[ the probabil- 
ity of success is asymptotically given by trs{nop(0)}. 
Therefore, except for those situations wherein the max- 
imum eigenvalue (in modulus) is unity, the most rele- 
vant condition to get a good efficiency is that the num- 
ber of required measurements is very low, which implies 
l-^i/Aol ^ 1, or, better, that the denominator in the 
threshold given in p5)) is high, i.e. In |Ao/Ai| » 1. There- 
fore, the peaks of ln|Ao/Ai| correspond to the maxima 
of the efficiency (i.e., the minima of the required num- 
ber of measurements). To better fix the idea, if we 



9 




10 



FIG. 2: The purity of the extracted state as a function of 
the temperature T (in units of ftSl/fca) and er. The other 
parameters are; 6 = 37r/4, x — 0, Q,/e = 10, 72/e = 0.1, 
71/72 = 0.95. 

ask that the target state is obtained with a precision 
Po — 0.99, since we have ln(4 — 1) + ln[po/(l — Po)] ~ 5.7 
we find that, in correspondence to those peaks wherein 
ln|Ao/Ai| ~ 4, the process requires one or two measure- 
ments when the system starts with an initial condition 
satisfying ||nojo(0)|| ^ R{p{Q)), which for instance is usu- 
ally the case for the maximally mixed state. 

In Fig. Efa), we consider the ideal case, while in 
Figs. [3Ub)-[3fd), we refer to non-ideal situations at zero, 
intermediate, and high temperature. The plots clearly 
show that the interaction with a nonzero temperature en- 
vironment negatively affects the efficiency, lowering the 
peaks and extending the valleys. Nevertheless, at zero 
temperature, various peaks are still present, and in fact, 
at zero temperature, the degradation with respect to the 
ideal case is not so dramatic. 



IV. SUMMARY 

Let us summarize the results reported in this paper. 
Putting a system in interaction with a repeatedly mea- 
sured one forces the former system onto a subspace, hence 
realizing, under suitable conditions, the extraction of 
pure states. In a more realistic situation, the two systems 
are interacting with their environment too, and therefore 
are subjected to dissipation. Such an interaction practi- 
cally reduces the chance to extract pure states. 

From the mathematical point of view, the main dif- 
ference between the two situations is represented by the 
fact that in the ideal case one extracts eigenvectors of a 
map onto a Hilbert space, whereas in the non-ideal case 
one extracts eigenvectors of a map onto a Liouville space. 



We have explored the general framework and studied a 
very simple physical system (two spins interacting with 
a bosonic bath) in order to bring to light fundamental 
features of repeated-measurement based extraction pro- 
cesses in the presence of dissipation. In Sec. IIII Bl we 
have shown that a mixed state is extracted instead of 
a pure state. Actually, this is what generally happens, 
especially at high temperatures. Nevertheless, with a 
zero-temperature bath, it is still possible to extract pure 
and almost pure states (see Fig. [TJb)) with still fairly 
good efficiency. Indeed, the efficiency, though negatively 
affected by the environment, is still good at zero temper- 
ature. 

Overall, we have considered the case wherein a very 
large number of measurements (evaluating the mathe- 
matical limit for an infinite number of measurements) is 
performed on the ancilla system, as clearly expressed by 
(fTU)) and pT|) . We conclude this paper expecting that in 
some cases a reduction of the number of measurements 
performed on the ancilla system entails an increase of 
the purity of the output state. See, for example. Fig. 21 
wherein we have plotted the purity of the resulting quan- 
tum state as a function of the number of measurements 
performed on the ancilla system, starting from the maxi- 
mally mixed initial state, with a particular parameter set. 
It is well visible that the purity, starting from its min- 
imum value (i), increases at the second measurement, 
and then decreases down to its asymptotic value. There- 
fore, in such a case, the highest value of purity is obtained 
for a smaller number of measurements {N = 2). We will 
discuss this aspect of our scheme in the next future. 
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FIG. 3: The quantity In |Ao/Ai| vs the parameters er and S/tt (x = for aU cases), in different situations: (a) in the ideal case, 
i.e. in the absence of interaction with the bath, and in the presence of interaction with the bath with (b) fceT/ftfi = 0.01, (c) 
kBT/hil = 1, and (d) ksT/hi} = 10. In aU cases, ratios between salient physical quantities are fixed as i}/e = 10, 72/e = 0.1, 
71/72 = 0.95. 
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N 

FIG. 4: The purity of the state of system S as a function of 
the number of measurements on X. The initial state is p(0) — 
0.5|t)s(T| + 0.5|i)s(i| and its purity is 0.5. Observe that the 
purity reaches the maximum value at A'^ = 2 and then decays 
to an asymptotic value. Here, the parameters are x = 0, 
9 = 2.25, and er = 7.82, and we have set kjiT/hO. = 10"^, 
n/e = 10, 72/e = 0.1, 71/72 = 0.95. 
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